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Abstract 



The two-fermion relativistic wave equations of Constraint Theory are re- 
duced, after expressing the components of the 4x4 matrix wave function 
in terms of one of the 2x2 components, to a single equation of the PauH- 
Schrodinger type, vahd for all sectors of quantum numbers. The potentials 
that are present belong to the general classes of scalar, pseudoscalar and 
vector interactions and are calculable in perturbation theory from Feynman 
diagrams. In the limit when one of the masses becomes infinite, the equation 
reduces to the two-component form of the one-particle Dirac equation with 
external static potentials. The Hamiltonian, to order 1/c^, reproduces most 
of the known theoretical results obtained by other methods. The gauge in- 
variance of the wave equation is checked, to that order, in the case of QED. 
The role of the cm. energy dependence of the relativistic interquark confining 
potential is emphasized and the structure of the Hamiltonian, to order !/(?, 

corresponding to confining scalar potentials, is displayed. 
PACS numbers : ll.lO.Qr, ll.lO.St, 12.40.Qq. 



2 



Typeset using REVT^iX 



I. INTRODUCTION 



The use of the manifestly covariant formahsm with constraints |l| in the two-body prob- 
lem leads to a Poincare invariant description of the dynamics of the system with the 
correct number of degrees of freedom. Furthermore, the potentials that appear in the cor- 
responding wave equations are calculable, in perturbation theory, in terms of the kernel of 
the Bethe-Salpeter equation, and therefore allow one to deal with quantum field theoretic 
problems. 

For two spin-0 particles, the system is described by two independent wave equations, 
which are generalizations of the individual Klein-Gordon equation of each particle, including 
now the mutual interaction potential. The compatibility condition of the two equations 
imposes certain restrictions on the structure of the potential and leads, in a covariant form, 
to an elimination of the relative energy variable. One then ends up with a final, manifestly 
covariant, three-dimensional eigenvalue equation that describes the relative motion of the 
two particles. This equation is very similar in form, except for kinematic energy dependent 
factors, to the Schrodinger (or Klein-Gordon) equation : it is a second order differential 
equation in the three spacelike coordinates (for local potentials) and therefore the usual 
techniques of nonrelativistic quantum mechanics are applicable to it. 

For two spin-1/2 particles, the system is described by two independent Dirac type equa- 
tions. Again, the compatibility condition imposes restrictions on the stucture of the poten- 
tials and eliminates the relative energy variable; but, because of the presence of the Dirac 
matrices, the reduction to a final eigenvalue equation is not straightforward. Except for cer- 
tain classes of interaction, as the pseudoscalar interaction and the vector interaction without 
temporal components, the reduction process is rather complicated and dependent on the way 
of eliminating the components of the spinor wave function in terms of one of them. Up to 
now, no single Pauli-Schrodinger type equation was obtained from this procedure, capable of 
describing general cases of interaction and all sectors of quantum numbers, although it was 
possible to obtain two coupled Pauli-Schrodinger type equations for the sectors of quantum 
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numbers with i = | j ± 1 1 and single equations for the sectors with i = j 0] . 

The purpose of the present paper is to show that it is possible, for general combinations 
of scalar, pseudoscalar and vector interactions, the latter being considered in arbitrary co- 
variant gauges, to reduce the wave equations describing two spin-1/2 particle systems to a 
single Pauli-Schrodinger type equation, valid for all sectors of quantum numbers. This is 
achieved by decomposing the 4x4 spinor wave function (of a fermion-antifermion system, 
say) along 2x2 components defined, in the cm. frame, by the basis of the matrices 1, 70, 75 
and 7o75, rather than the usual basis (1 ± 70) and (1 ± 70)75. 

The above reduction is obtained in the case of quasilocal approximations of the interac- 
tion potentials, i.e., for potentials that are functions of the relative coordinates and, eventu- 
ally, of the cm. energy. This does not cover the more general case of integral operators (in 
the three spacelike coordinates), but it is generally admitted that a local expression of the 
potential may provide a valid zeroth order approximation of the total interaction, from which 
one may develop perturbative calculations for the nonlocal effects. Furthermore, an appro- 
priate cm. energy dependence of the potential takes into account the leading contribution 
of the nonlocal effects. 

The final Pauli-Schrodinger type equation that is obtained provides explicit and simple 
means of controlling the relativistic effects of each type of interaction. 

The contents of the paper are summarized as follows. In Sec. II, we review the main 
aspects of the relativistic wave equations of Constraint Theory for a fermion-antifermion 
system. The calculation of the potentials from one-particle exchange Feynman diagrams is 
presented in Appendix A. 

In Sec. Ill, we reduce these wave equations to a single two-body Pauli-Schrodinger 
type equation, for the case of a general combination of scalar, pseudoscalar and vector 
interactions. Some technical details of the reduction process and useful formulas for the 7- 
matrices and the spin and orbital angular momentum operators are presented in Appendix 
B. 

In Sec. IV, we study the limit of the above equation when one of the particles becomes 
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infinitely massive. It reduces to the Paufi-Schrodinger form of the one-particle Dirac equation 
in the presence of the static potential created by the heavy particle. In Appendix C, we 
show that, for the class of ladder and crossed ladder type diagrams, it is only the one-particle 
exchange diagram which contributes in this limit to the expression of the potential. 

Section V deals with the nonrelativistic limit, to order of the above two-body Pauh- 
Schrodinger type equation. Most of the known theoretical results obtained by other methods 
are reproduced. In the case of QED, its gauge invariance, to order 1/c^, is checked by 
considering the photon propagator in arbitrary covariant gauges. The role of the cm. energy 
dependence of the relativistic potentials is emphasized and the structure of the Hamiltonian, 
for confining scalar potentials, ready for comparisons with quantum field theoretic results, 
is displayed. 

Summary and concluding remarks follow in Sec. VI. 
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II. THE TWO-BODY RELATIVISTIC WAVE EQUATIONS 
OF CONSTRAINT THEORY 



For a system of two spin-1/2 particles, composed of one fermion of mass mi and one 
antifermion of mass m2, say, Constraint Theory imposes two independent wave equations, 
which are generahzations of the Dirac equation satisfied by each particle in the free case. 
These wave equations have not a unique form and can be modified by wave function trans- 
formations, but a convenient form, where several properties can be easily read off, is the 
following one : 

( 7i.pi - mi) ^ = ( - 72.P2 + ma) r ^ , (2.1a) 
(-72-^2-^2)^ = ( 7i.pi + mi) \/ § . (2.1b) 

Here, \l/ is a sixteen-component spinor wave function of rank two and is represented as a 
4x4 matrix : 

^ = ^a,,a,{Xl,X2) (c^i, ^2 = 1, • • • , 4) , (2.2) 

where 01(02) refers to the spinor index of particle 1(2). 71 is the Dirac matrix 7 acting in 
the subspace of the spinor of particle 1 (index ai); it acts on \^ from the left. 72 is the Dirac 
matrix acting in the subspace of the spinor of particle 2 (index 0:2); it acts on \E' from the 
right; this is also the case of products of 72 matrices, which act on \l/ from the right in the 
reverse order : 

7lM^ = (7M)ai/3i^/3ia2 ' 72/.^ = ^ai/Jj (7/.)/32a2 , 
72/.72;.^ = ^a^|32ilul^,)f32a2 , = ^baa^^a/s] (o = 1, 2) . (2.3) 

In Eqs. (2.1) pi and p2 represent the momentum operators of particles 1 and 2, respectively, 
is a Poincare invariant potential. 

Equations (2.1) must be compatible among themselves. This is enforced in two steps. 
First, one multiplies Eq. (2.1a) by (7i.j9i + mi) and Eq. (2.1b) by (—72.^2 + ^2) and 
subtracts the two equations from each other. This yields the constraint : 
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[{pl-pl)-{ml-ml)]^ = 0, (2.4) 

which allows one to eliminate the relative energy variable in a covariant form. For eigen- 
functions of the total momentum operator P, the solution of Eq. (2.4) is : 

^ = e-iP-^ e-^(^i - ml)Rx/{2P') ^^^t^ ^2.5) 

where we have used notations from the following definitions : 

-P = Pi+P2, P=^(Pi-P2), M = mi + m2, 

X = i(x,+..). . = x,-x.. (2.6) 

We also define transverse and longitudinal components of four-vectors with respect to the 
total momentum P : 

^ = <li-iQ-P)P,, P.^Pjv^, 

qL = q.P : Pl = Vp^ . (2.7) 

This decomposition is manifestly covariant. In the cm. frame the transverse components 
reduce to the three spacelike components, while the longitudinal component reduces to the 
timelike component of the corresponding four-vector. (Note that x^^^ = — in the cm. 
frame. ) 

Second, one multiplies Eq. (2.1a) by (— 72.p2 — '^'2) and Eq. (2.1b) by (71 .pi — mi) and 
subtracts the two equations from each other. Taking into account the constraint (2.4), one 
ends up with the new constraint : 

[pI-pI ,V^] ^ = 0, (2.8) 

which concerns the potential V. By noticing that pi—pi — '^PlPl, this equation is satisfied 
by demanding that V, which is Poincare invariant, be independent of the longitudinal relative 
coordinate xl 

V = V^(x^, Pl, p^, 71, 72) ■ (2.9) 



No other constraints are found. Equations (2.5) and (2.9) show that the internal dynam- 
ics of the system is three-dimensional, besides the spin degrees of freedom, described by the 
three-dimensional transverse coordinate . 

To find the connection of Eqs. (2.1) with the Bethe-Salpeter equation one projects 
the latter on the constraint hypersurface (2.4) by expanding the right-hand side of the 
equation on that hypersurface [0]. One ends up precisely with Eqs. (2.1), where now the 
potential V is calculable, in perturbation theory, in terms of Feynman diagrams. 

The relationship between the potential V and Feynman diagrams can be summarized 
by the following Lippmann-Schwinger-Quasipotential type ||8|-[15| equation (momentum in- 
tegrations, as well as the total momentum P are not explicitly written) : 

V - f - VGof = , (2.10) 



T{p,p') 



C{p),C(p') 



2Pl 
where : 

i) T is the off-mass shell fermion-antifermion scattering amplitude (amputated four-point 
connected Green's function); 

ii) C is the constraint (2.4) : 

C{p) = (pI-pI) - {ml - ml) = 2PlPl - {ml - m^) ^ ; (2.11) 

in Eq. (2.10) the external momenta of the amplitude T are submitted to the constraint C; 

iii) Go is defined as : 

^0(^1,^2) = SM H-V2) Ho , (2.12) 

where Si and 5*2 are the propagators of the two fermions, respectively, in the presence of the 
constraint (2.11), and Hq is the Klein-Gordon operator, also in the presence of the constraint 
(2.11) : 



Ho = {Pl - m\] 



= {pl-mD 

c 



^ = ^-lirnl + ml)+^^^^^L^+p^-. (2.13) 



The series expansion of Eq. (2.10) involves not only the usual Feynman diagrams of the 
amplitude T, but also additional diagrams, which we call "constraint diagrams", coming 
from the iteration of the VGq term; they are obtained from the reducible diagrams by the 
replacement in a box diagram of the fermion and antifermion propagators by the factor 
Go (2-12) together with the (5-function of the constraint (2.11). It is the presence of these 
diagrams, as well as of the cm. energy factor in the function S{C), that cures the Bethe- 
Salpeter equation of some of its deficiencies (See also Appendix C.) [We have included 
the cm. energy factor 1/(2Pl) in the definition of T, and therefore, by imposing on Go 
the constraint (2.11), the integrations relative to Go in Eq. (2.10) may be considered as 
three-dimensional.] 

In general V is an integral operator in x^. However, in one-particle exchange diagrams 
V turns out to be a function of x'^ and Pl, with definite dependences on the 7-matrices. 
It would also be meaningful to approximate, as a zeroth order approximation, multiparticle 
exchange contributions by appropriate local functions and thus to use for V quasilocal 
expressions (in and Pl). Furthermore, in perturbation theory, the cm. energy (Pl) 
dependence of the potential takes into account the leading contribution of the nonlocal 
effects. Therefore, we shall confine ourselves in the present paper to quasilocal expressions 
of V. The explicit expressions of V in lowest order of perturbation theory, for the scalar, 
pseudoscalar and vector interactions, are presented in Appendix A. 

The scalar product associated with the solutions of Eqs. (2.1) can be constructed by 
means of tensor currents of rank two, satisfying two independent conservation laws, with 
respect to xi and X2 0. Another way of proceeding consists in using the integral equation 
of the Green's function [jl2|,|15[. From either method, one finds for the norm of ip [Eq. (2.5)] 
the formula (in the cm. frame) : 



(The potentials V that are considered in the sequel commute with the matrix product 




(2.14) 



71L72L-) 
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This formula shows that for energy independent potentials (in the cm. frame) the norm 
of ip is not positive definite. In order to reach such a situation, it is sufficient that the 
potential V satisfy the inequality 

^Tr{V^V) < 1 , (2.15) 

which, in turn imposes restrictions on the high order (multiparticle exchange) contribu- 
tions to V. Once condition (2.15) is satisfied, one is allowed to make the wave function 
t r ansf ormation 

§ = [i-v^v]'^' ^ (2.16) 

and to reach a representation of the wave equations (2.1), where the norm for energy inde- 
pendent potentials is the free norm. 



In this connection, it was observed by Crater and Van Alstine [|T6[ that the parametriza- 
tion 

V = tanhV (2.17) 

satisfies condition (2.15) and also has the property that in the new representation (2.17) one 
obtains familiar Dirac type wave equations, where each particle appears as placed in the 
external potential created by the other particle, the latter potential having the same tensor 
nature as potential V of Eq. (2.17). 

We shall henceforth adopt parametrization (2.17) for potential V and, according to Eq. 
(2.16), shall introduce the wave function transformation : 

^ = (coshy) ^ , (2.18) 

V being assumed to be hermitian. The norm of the new wave function \1' then becomes (in 
the cm. frame) : 

J d^x Trjv^t [1 + 4^^Q^2oPo ^] = 2Po • (2.19) 
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(The relationship between ^ and t/j is the same as in Eq. (2.5).) 

The new potential V will be chosen in the major part of this paper as a general combi- 
nation of scalar, pseudoscalar and vector potentials : 

T T 

V = V, + 715725^^3 + 7^72 ( qO^V^ + gluU, + ) , (2.20) 

where the potentials V^, V2, etc., are functions of x^"^ and P^. (The index notations of 
the potentials will be clarified in Sec. III.) The dependence of the vector potentials on the 
gauge chosen for the photon propagator in lowest order of perturbation theory is displayed in 
Appendix A. Also notice that the difference between potentials V and V shows up starting 
from third order diagrams. 

When parametrization (2.17) is adopted and the change of function (2.18) is used, Eqs. 
(2.1) take the form : 

( 7i-Pi — ^1) coshy %l) — (— 72.p2 + ^77.2) sinhy -0 , (2.21a) 
(— 72-P2 — "^2) coshy -0 = (7i-Pi + ?«i) sinhV^-^, (2.21b) 

where, according to Eq. (2.5), ip{x^) is the internal part of the wave function ^. Also, 
because of the constraint (2.4) (or (2.11)), the longitudinal components of the momentum 
operators are replaced by their eigenvalues determined in terms of Pl [Eq. (2.7)] and the 
masses : 

p,, = ^+(!li=^. (2.22a) 
p., = ^-M^. (2.22b) 

By using the expression of V [Eq. (2.20)] and the properties of the 7-matrices (see 
Appendix B) , one can bring the functions cosh V and sinh V to the left of the Dirac operators 
in Eqs. (2.21). If one defines V = 2Vi — V, then, upon multiplying Eq. (2.21a) by coshV^ 
and Eq. (2.21b) by — sinh V and subtracting the two equations from each other, one obtains 
a Dirac type equation for particle 1 as if it was placed in the presence of external potentials 
created by particle 2 : 
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2 2Pr ^ '^^ 2 2M 



+ e" 



-2C/4 



I V' = 0, 



T T T , ,T 

- 2ihe-'^'^^^lx^ + 7il72lT>2 + 7i5725l>3 + 7^ ■72'^f>4 + J,' t 



(2.23) 

where we have defined : 

dF 

F^S^.- (2.24) 

The Dirac type equation for particle 2 can be obtained from Eq. (2.23) by the replace- 
ments : pi — p2, X — >• X, mi 7712, 7i ^ 72 ■ 

Equation (2.23) shows that the scalar potential acts essentially on the mass terms, the 
timelike vector potential acts on the energy terms, the spacelike vector potentials act on the 
transverse momentum dependent terms, while the pseudoscalar potential acts only through 
spin dependent terms. These results are consistent with the interpretation of Eq. (2.23) 
from the external potential standpoint. In particular, we did not find, in passing from Eqs. 
(2.21) to Eq. (2.23) new types of interaction. This is not a trivial property. Had we used 
a different parametrization than that provided by Eq. (2.17), for instance by decomposing 
V, rather than V, according to Eq. (2.20), we would have found that the original vector 
potentials induced in the analogue of Eq. (2.23) additional axial vector terms. 
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III. REDUCTION OF THE WAVE EQUATIONS 
TO THE PAULI-SCHRODINGER FORM 



In order to bring Eqs. (2.21) into a more transparent form, we decompose the 4x4 
matrix wave function on the basis of the matrices 1, 7^,, 75 and 7L75 by defining 2x2 
matrix components : 

4 

i=l 

Equations (2.21) are then projected in the above subspaces. We shall first consider 
the case of more general potentials than those introduced in Eqs. (2.20) and develop our 
method of calculation for this general case. They have the property of being functions of 
products of 7i and 72 matrices in equal number. This feature is obvious for one particle 
exchange diagrams (for parity conserving interactions) , where the matrices 71 and 72 couple 
symmetrically to the propagator of the exchanged particle. However, when general vertex 
corrections are considered, one finds expressions in which some terms may not have the same 
number of 71 and 72 matrices, although globally the total expressions are symmetric in the 
exchanges 1 2. Such terms will not be considered here. Then, the most general (parity 
and time reversal invariant) potential we may consider has the following decomposition on 
the basis (3.1) : 

V ^ Vi+ 71L72L V2 + 715725 ^3 + 7il71572L725 V4 , (3.2) 

where the potentials Vi {i = 1, . . . , 4) may still have spin dependences. 

Spin operators, which act in the 2x2 component subspaces, are defined by means of the 
Pauli-Lubanski operators : 

Wis = Wis = ' Ws = W^s + . (3.3) 

They also satisfy the relations : 
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hP hP 

'flLWiSa = ^7^715 , l2LW2Sa = ^7^725 ■ (3.4) 



(Note that the Ws are transverse vectors.) 
We introduce the operators : 
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w 



2 1^,^.1^25 ^ -7281-82 , (3.5a) 



2 2 WiS-x'^W2S-x'^ cm. 4 (Si.x)(s2.x) 



In the preceding relations we also displayed the corresponding cm. expressions in terms of 
the usual spin operators. 

Then, the potentials Vi of Eq. (3.2) can be decomposed as : 

Vi ^ Vi + wUi + wuTi (i = l,...,4) , (3.6) 

where Vi, Ui and Tj are functions of and Pl- 

For the combination of potentials we introduced in Eq. (2.20), we have to reexpress the 
transverse matrices 7^ according to Eqs. (3.4). By comparing Eq. (2.20) with Eqs. (3.2) 
and (3.6) we identify the indices used in Eq. (2.20). 

In order to project Eqs. (2.21) onto the subspaces of the components ■^j [Eq. (3.1)], it is 
useful to introduce projection matrices for these subspaces. These are : 

Vi = ^(1 + 7il72l) (1 + 715725) , ^2 = ^(1 + 7il72l) (1 -715725) , 

1 1 

V3 = ^(1 - 7il72l) (1 + 715725) , = ^(1 -7il72l) (1 -715725) ■ 

(3.7) 

They satisfy the relations : 

ViVj ^ s.jVj , v^^, = 6^,r, (i,j = i,...,4) . (3.8) 

(The r's are defined in Eq. (3.1).) 

The general potential V [Eq. (3.2)] can be expressed in terms of the projection matrices 

V : 
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i=l 



Qi ^ Ai + wBi + wi2Ci ; (3.9) 

the relations between the sets A^, B^, Ci and V^, C/j, Tj [Eq. (3.6)] are easily established. 
For the potential (2.20) we shall consider below, we have the identifications : 

ai = 14 + + "Vs + wUi + ^12^4 , 

a2 ^ Vi + V2 - - WU4 - W12T4 , 

as ^ Vi — V2 + V3 — WU4 — W12T4 , 

04 = 1/1 - 1^2 - ^3 + WU4 + W12T4 . (3.10) 

The projectors (3.7) are mainly needed for the evaluation of the exponential functions 
of V, present in Eqs. (2.21). One has the simple property : 

4 4 

exp(5]a,n) = E^^e«^ ■ (3-11) 

1=1 i=i 

Equations (2.21) then yield the eight coupled (compatible) equations for the four com- 
ponents ipi {i — 1, . . . ,A) : 

e-«i ^1 - M e«2 ^2 + ^Ws.p e'^^ = , 
Pl nPL 

(jlizJI^ _ M e«i V'l - j^^{W,s - W2s).p e'^^ V'4 = , 

e«3 ^3 - ^^^^ e-^ ^4 + ^Ws.p e-^ = , 
-Pl M KPl 

e«4 - K_l!^ e-«3 V.3 - ^(ly,^ - W2s).p ^^2 = , 
Pl ^1 - M— !M e^«2 + ^{W,s - W2s).p e"^ ^3 = 0, 

IVl rirL 



Pl ^2 - e-«i V'l - ^^s-P e^^ ^4 = 0. 



(3.12) 
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By using the property {Ws.p){Wis — 1^25).^ = 0, one easily deduces from Eqs. (3.12) 
the following four equations : 

Ws.p {Pl e~^3 ^3) ^ Ws.p (M e«4 ^ 
(Wis - W2s).p {Pl e-"4 = ^Wis - W2s).p (M e«3 ^3) , 
{W,s - W2s).p {Pl ^^;,) = {W,s - W2s).p i^^^^Z^ ^-^^ , 
Ws.p {Pl V'l) = Ws-P ( ^^'^^^^ e-«2 , 

(3.13) 

which are useful during intermediate stages of calculation and for consistency checks. 

Among the four components it is the combination ■?/'3 H-'?/'4 which is the dominant one 
in the nonrelativistic limit. The quantum numbers of the state are therefore determined by 
■03 and ■04. For this reason one must eliminate from Eqs. (3.12) ■0i and ■02 in terms of ■j/'a and 
■04 and finally keep ^03 or ■04. We shall keep the component ^03 and express the eigenvalue 
equation in terms of it; the calculations with -04 parallel those with 03 and we shall indicate 
at the end the substitutions to use to obtain the eigenvalue equation with ^04. 

The third and seventh of Eqs. (3.12) allow us to express 0i and 02 in terms of 03 : 

^ g(ai - 2/i) _|_ I _ g-2ai ^^^^ _ ^^^^ ^ ^03 

+ ^"^^"^ + "^S.P e-«^ 03 } , (3.14a) 

= e(«2 - 2/i) I + e-2a2 ^^^.^ ^-03 ^3 



mr — m 



p2 

where we have defined : 



I) g-2(ai + 02) (^^^ _ ga3 ^3 I ^ (3_i4b) 



= 1 - l!li^^e-2(«i + «2) (3.15a) 



^ _ (mj - m^)^ 4(v, + 1^2) 
M2P2 

the second equation, (3.15b), being valid only for the case of potential (2.20). 
One then obtains two independent equations for -03 and -04 : 
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(3.15b) 



MPl e~^4 




(3.16a) 




Finally, eliminating ■04 through Eqs. (3.16) one obtains the eigenvalue equation with ■03: 



Had we expressed ■^i and 1/^2 in terms of 1/J4 and eliminated ij^s in terms of ^04, we would 
have obtained the eigenvalue equation with ^04, which can be read off Eq. (3.17) with the 
following substitutions in the terms containing the spin operators : —04 and 



Equation (3.17) is a second order differential equation and has the form of a generalized 
Pauli-Schrodinger equation. It is possible to bring it into a more standard form by moving 
the momentum operators to the utmost right and having recourse to wave function transfor- 
mations. We shall undertake these operations for the case when potential V takes the form 
(2.20) (or (3.10)) (scalar, pseudoscalar and vector potentials), the calculations, as well as the 
final expressions, being relatively simple. For simplicity of notation and a more transparent 
reading of the formulas, we shall express the resulting equation in the cm. frame. We use 
the notations : 



p2g-(a3 + 04) _ ^2g(a3 + 04) 

-(-^)2 e«4 Ws.p e-^^ { e"2«2 ^ g-oa 



L 




(3.17) 




|)/P2 ^ (ml - ml)/M\ 
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S = si + S2, L = xxp, J = L + S, Su = . (3.18) 

ft X 

Upon bringing, in Eq. (3.17) and for the case of potential (2.20) (or (3.10)), the momen- 
ntum operators to the utmost right, one finds that the quadratic terms in p sum up to yield, 
apart from a multiplicative potential dependent factor, the operator p^. However, one also 
meets terms proportional to i/ix.p. To get rid of the latter terms, it is necessary to use the 
change of function : 

= e(2^4 + 2T4 + /i) I 1 (1 ^ g^^^ ^2V, ^ _ ^2V2 I _ (319) 

After some lengthy calculations (see Appendix B for useful formulas that are utilized for 
the spin and orbital angular momentum operators) one ends up with the following eigenvalue 
equation : 
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P^AV2 _ MljVi _ {ml-mlf , (^i " ^If -4V2 

4 4 4M2 ^ 

- - 4h^^^ {v; + v^ + v^-ui + Ti + h')^ 

+ eh" ( v; + \/2' + v^3 - ^1 + + /i' ) + Ah^^^ ( v;" + + v^g" - + + h" ) 



48^ ( Vi + Vi + Vi + U'^ + T'^ + h' ) {1 + A^^U'^ ) 

-(C/i + 2x=^C/^') 
(S.x)2 



' (e^^^ - 1) 



4x2 



+ 4^^:^=^ ( V; + V^ + 2U'4 + Ti + h' ) {1 + 4:k.\V^ + U^) ) 

2x^(^3" + 

(S.x)2 



+ 4( L.S + S' 



2x2 

2T4 



sinh^(yi - ^2) 



^ 'e^^^ _ 1) cosh 2(^1 -^2) 



2x2 



+ {Vl + + 2^/4 + h' ) cosh 2(^1 - V2) - (V; - V^) sinh2(Vi - V2) 



+ 4[|^l!,L.Sle2r4 

h x^ 



^ ^e2^4 _2) sinh2(yi -^2) 



2x2 



- {Vl + V2 + 2U'^ + h' ) sinh 2(1^1 - V2) + (l^/ - 1^2) cosh 2{Vi - V2) 
— 4:{ml — ml) 



1 


1 






1 


1 


P2 





. 



(3.20) 



We recall that the functions Vi, V2, V3, U4 and T4 correspond to the scalar (Vi), timelike 
vector {V2), pseudoscalar (^3) and spacelike vector (C/4 and T4) potentials. The function h 
is defined in Eq. (3.15). We have also defined : 

dF „„ d^F 



F' = 



F" = 



(3.21) 



9x2 ' ^ - (5x2)2 

Equation (3.20) is the main result of the present paper. It corresponds to the reduction 
of the initial covariant wave equations (2.1), by means of the parametrizations (2.17) and 
(2.20) and of the wave function transformations and decomposition (2.18), (3.1), (3.19), to 
the Pauli-Schrodinger form, in which ^3 is a 2 x 2 matrix wave function. 
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The above reduction can be realized with respect to the component ^4 as well. In this 
case the wave function transformation (3.19) is replaced by the following one : 

^4 = e(2t^4 + 2T4 + /i) I 1 ^ ^2V2 ^ hi _ 3^2) e^^i } 04 . (3.22) 

The equation satisfied by 04 is the same as Eq. (3.20), except a global change of 
sign in front of the terms proportional to the commutators [ (S.x)^/(^^x^) , L.S ] and 
[ (S.x) V(^V) , L.(si - S2) ] . 

In order to reconstitute the whole wave function %l) from V's, it is also necessary to calculate 
-04 from ■^s (or conversely, V's from ■^4 if the latter is first obtained). To this end, the first two 
of Eqs. (3.13) kinematically yield the desired relations in the sectors of states corresponding 
to the quantum numbers j — £ . However, in the sectors with quantum numbers j ^ i 
these equations do not allow one to obtain kinematically the relationship between ip^ and 
■04 . It is therefore necessary to go back to the two coupled equations (3.16). Either of these 
equations can be used to express 1^4 in terms of 1^3 . As they stand, these equations involve 
second order differential operators. In order to obtain a relationship involving at most first 
order differential operators, one eliminates the operator (see for this Eq. (B8)) by the 
use of Eq. (3.17). The relation giving ^04 in terms of V's is presented in Appendix B [Eq. 
(B9)] where we also explain how to express ■ip3 in terms of ■ip4^ . 

For completeness, we rewrite, in Appendix B [Eqs. (BIO)], Eqs. (3.14) giving ■0i and -02 
in terms of ■03, in more explicit forms when potential V has the form (2.20). 

Equation (3.20) can be solved by decomposing the wave function 03 along radial wave 
functions and the generalized spherical harmonics, yjl^ , eigenf unctions of the operators 
J"^) Jz) L"^, S"^ . 

When all potentials are present, and in the unequal mass case, there are only two sectors 
of quantum numbers where 03 is an eigenfunction of the latter operators : i) s = 0, i — 
0, j — ] u) s — 1, i — 1, j — . For the other sectors 03 is a combination of different 
eigenfunctions : s = and s = 1 for the sectors with i = j ; i = \j ± 1\ for the sectors with 
s = 1 . In the equal mass case, a simplification occurs, in that for the sectors with i — j, 
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the subsectors with s = and s = 1 disantangle. A similar simphfication also occurs in the 
unequal mass case when Vi and V2 are absent, or when V1 + V2 = (or constant). Complete 
diagonalization occurs only for harmonic oscillators. 

The presence in Eq. (3.20) of the commutators [ (S.^)^ /{h^x^) , L.S ] and 
[ (S.x)^/(/i^x^) , L.(si — S2) ], might suggest that some of the energy eigenvalues are not 
real. These operators connect only off-diagonal matrix elements of the sectors i = |j ± 1| . 
As is known, such matrix elements contribute only with the modulus squared of their values, 
and therefore the resulting energy eigenvalues are real. This also explains why, on solving 
either of Eq. (3.20) or the equivalent equation of 04 , where the signs in front of the above 
commutators are changed, one obtains the same energy eigenvalues. 

In this connection, it is worth emphasizing that if, instead of considering the complete 
norm (2.19), where in any event the hermiticity property of the energy operator is manifest 
provided V is superficially hermitian, one tries to construct the scalar product or the norm 
for 03 directly from Eq. (3.20), one has to consider as the adjoint of 03 the function rather 
than 03 , precisely because of the antihermitian nature of the abovementioned commutators. 
If Eq. (3.20) is written in compact form as ^303 = , then following standard methods 
P JT^JT3| , one ends up with a norm proportional to / (i^x Tr |p| 03 , the hermitian 
conjugate of which is nothing but / d^yi Tr 03 |p| 04 , because L3 = L4 , where L4 is the 
wave equation operator of 04 , satisfying 1^404 = . 

Finally, we observe that when the potentials Vi and V2 are absent, then 03 and 04 
satisfy the same eigenvalue equation, and yet they are not identical in general (see Eq. 
(B9)). This means that the energy spectrum has then degeneracies coming from the sectors 
with i = |j ± 1| . These degeneracies increase in the case of potentials yielding harmonic 
oscillators. 
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IV. INFINITE MASS LIMIT 



In this section we shall study the limit of Eq. (3.20) when one of the particles, 2, say, 
becomes infinitely massive. It is known that the Bethe-Salpeter equation does not have a 
simple behavior in this limit. The origin of this difficulty lies in the fact that the kernel of the 
Bethe-Salpeter equation contains only irreducible (crossed ladder) diagrams, which do not 
have suitable limits. In order to reproduce in this limit, from the Bethe-Salpeter equation, 
the one-particle Dirac equation in the presence of the potential created by the infinitely 
massive particle, one must take into account contributions coming from the infinite number 
of crossed ladder diagrams [|l7i[T8| . 

The situation is much simpler for the wave equations of Constraint Theory. Here, the 
potential is calculated from the diagrams of the off-mass shell scattering amplitude, to which 
one adds "constraint diagrams" , obtained from reducible diagrams with the use of constraint 
(2.4) [Eq. (2.10)]. To a given formal order of perturbation theory, the sum of all ladder 
and crossed ladder type Feynman diagrams (without vertex corrections) does have a simple 
limit : the sum of the n^^ order diagrams (in coupling constant squared) contains the factor 



S{kio) ■ ■ ■ 6{kn~io), where the fc's are the internal momenta of the loops {n > 2) ||T7| , p!8[| . This 
term is then cancelled by the dominant term of the "constraint diagrams" of n^^ order ||^ 
(see Appendix C). Therefore, in the limit m2 — oo it is only the one-particle exchange 
diagram (including vacuum polarization) which survives and leads to the one-particle Dirac 
equation with the corresponding static potential. 

The expressions of the potentials calculated from such diagrams are presented in Ap- 
pendix A. They essentially have the energy dependence of the form 

^ = ^ V^'^ , (4.1) 
and, therefore, when m2 goes to infinity their behavior is : 



lim V = +0(^). (4.2) 

m2^oo 2m2 mi 



22 



The contributions coming from higher order diagrams are of order 0{l/ml). If one 
replaces in Eq. (2.23) the various potentials by their leading behaviors (4.2), one immediately 
finds that Eq. (2.23) reduces to the Dirac equation of particle 1 in the presence of the external 
static scalar and timehke vector potentials v}^^ and V2^\ It is, however, also instructive 
to study the above limit in Eq. (3.20), where the contributions of each potential are more 
involved. 

To this end we expand each potential Va {a — 1, . . . ,A, V4 representing here either U4 or 
T4) up to its next to the leading term: 



The total energy Pl is decomposed as: 

Pl ^ m2 + E . (4.4) 
The terms inside the brackets containing Vi and V2 have the following hmit: 



lim 

m2— >oo 



El JV2 _ Ml AVi _ K - ^i)' -Wi , (^1 - miy 

4 4 ^ 4P^ 



(4.5) 



Notice that the terms W^^^ have disappeared, although the brackets were containing terms 
of the order of m|. It is also evident that the factor exp (4(C/4 + ^4)) in front of the brackets 
in Eq. (3.20) will contribute as unity. 

In the remaining terms of Eq. (3.20), besides the kinetic energy term, it is only those 
containing h', h" and exp(— 2/i) which will survive. They have the limits: 

m2^oo rn2 1 / J ' 

-(0)/ , T^(0) 



^2-00 2[{E + V^^^) + (mi + V^"') 



After averaging over the spin states of particle 2, one ends up with the following equation 
(where (p designates the limit of (p^/ 1712) ■ 
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- Ah^x^h'^ + m^h' + AU^^^h" + 8/i'L.si = . (4.7) 

This equation is nothing but the reduced form of the Dirac equation in the presence 
of the static scalar and timehke vector potentials after the decomposition — (^~) 
the change of function '?/'_|_ = {E + ¥2^'' + mi + arc used (in the basis where 70 

is diagonal). The latter is the limit (up to a multiplicative constant factor proportional to 
7712 ^^^)) when m2 goes to infinity, of the change of function (3.19). 
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V. NONRELATIVISTIC LIMIT TO ORDER l/C"^ 



We study, in this section, the nonrelativistic hmit of Eq. (3.20) to order 1/c^, the gauge 
transformation property of the resulting Hamiltonian for the electromagnetic interaction 
case and the structure of confining scalar potentials. 

A. General case in perturbation theory 

To recover the explicit c dependence of the various terms in Eq. (3.20) (where c had been 
set equal to 1) we must adjoin a c factor to each of the terms fi, p, Si, S2, L and a factor 
to mi and 777,2. On the other hand, at leading order, the potentials have, in perturbation 
theory, the dependence on of the type 1/{2Pl) (see Appendix A). Therefore, we expand 
the potentials as follows : 

K = ^ V}'^ + ^ V^'^ + 0{PE') (a = 1, . . . , 4) . (5.1) 

Here, V^*^°Y(2-Pl) is the contribution of the one-particle exchange diagram (or an effective 
representation of it) and V^^y{2Pl) is the leading contribution coming from the higher 
order diagrams. (In the perturbation theory of the present formalism, it is the two-particle 
exchange diagrams which provide the latter term.) 

To carry out the expansion, we split the total energy into the total mass and the binding 
energy E : 

Pl ^ M + E . (5.2) 

The term E, which is an eigenvalue, is brought to the utmost right and replaced there 
by the Hamiltonian H. In expressions containing the term H^, which enters only in 0{c~^) 
quantities, one replaces H by its nonrelativistic expression. 
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The resulting Hamiltonian has the the following expression (in the cm. frame) : 
^2 11 1 



H 



P_ + 1/(0) _ 1/(0) _ t ( ^ + fD2x2 



2mim2 



[p^ uf^+n 



<0) 



{ml 



1 



yiO) _ y(U) y2 ^ 



2M 

4mim2 



(0) n2 



M 



[p^ 



(0) 



+ 



(0) 



+ 



mim2 

1 fS.x)2 



mim2 



2mim2 



mim2 
1 

mim.2 

1 . 1 
2 



L.S 



(0)/ 



tP (mi - m2) 



+ 



2mim2 



(y/°^' + y2'"^0 



(0)/N 



2 

TTli 



ml 



(5.3) 



Here, /i = mim2/M, AV is the Laplacian of V, AV = 6V' + A-x^V", and the derivatives 
are with respect to . [ , ]+ is the anticommutator. To obtain Eq. (5.3) we also used the 
wave function transformation : 



03 



1 



1 



1 - ^rr/^o) + t(o)i - 



y{o) _ {mj + m^ - mim2) ^^(p) 



2mim2M 



(5.4) 



to bring the Hamiltonian into a manifestly hermitian form. The term containing 
[ (S.x)^/(^^x^) , L.S ] in Eq. (3.20), which is formally of order c~^, does not actually 

contribute, for it vanishes in diagonal matrix elements of the eigenvalues of the nonrelativis- 
tic Hamiltonian. 
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B. Electromagnetic interaction 



We now consider the case of the electromagnetic interaction. It is known that bound 
state energies are gauge invariant |19,2^. We check this property to order c~^. We set in 
Eq. (5.3) Vi = V3 = 0. Furthermore, considering the photon propagator [Eqs. (A6) and 
(A7)] we deduce that f/4 must be spht into two terms, V2 and a pure gauge term U; T4 is 
then equal to 2y^U' : 



V2 + U 



jj, _ 



(5.5) 



After this splitting is introduced, the Hamiltonian (5.3) for the electromagnetic interac- 
tion takes the form : 



H 



2/i 



1 1 

+ 



m2 



1 



( v^'^ + f7(°) ) p2 



+ 
+ 

+ 
+ 



17111712 

2ih 

17111712 

4:17111712 

2mim2 
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2LA(o)' 



2\2 



11111112 



X^XjPj^Pj 



^(0)' + 2f7(°)')x.p + 



ih 



17111712 



x.p 



A 



{6l7li17l2 



m: 



2i7iim2 
1 



2^ v^'^ + 3f7(°) 



2M 



( -Uvi'^ + V}'^' ) - 

171x1712 2 

nil ^2 77117712 



1 

M 

1 (S.x)2 



(0), 



mim2 ^ 

- ^ ( A - — ) L-(si 

2 mi 1712 



1 



S2) V, 



(0)/ 



(5.6) 



We notice that the gauge function [7^°^ has disappeared from the spin dependent terms. 
To prove gauge invariance of the energy corresponding to the Hamiltonian (5.6), we must 
show that a unitary wave function transformation removes the arbitrary gauge function {7*^°^ 
from H and brings it into its form of the Feynman gauge (in which U^^^ = 0), which is taken 
here as the reference gauge. 

It can be checked that the following wave function transformation 
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^ ( 1 + iG ) , 

G = J- ( f/W p.x + x.p t/W ) ^ ^ [ , r' [7(°) rfx^ ] , (5.7) 

where (pp is the wave function in the Feynman gauge, removes from Eq. (5.6) the existing 
U^^^ dependent terms. [Actually, it can be shown, by using the connection of the Constraint 
Theory wave function with the Bethe-Salpeter wave function 0, that transformation (5.7) 
is the three-dimensional reduction, to the leading order of the present approximation, of the 
gauge transformation of the Bethe-Salpeter wave-function 

The transformation (5.7) makes appear, however, in Eq. (5.6), a new U^^'^ dependent 
term, coming from the commutation of G with V2^^ and equal to 2x.'^U''^''V2^^'/M. This term 
must be compensated by the gauge dependence of the second order term V2^\ present in 
Eq. (5.6). If we designate by the expression of ¥2^^ in the Feynman gauge, then ¥2^^ is 
related to V2F by the relation: 

^(1) = V^P + 2x2 ^(0) ^(0), _ (5_g^ 

The expression of V2P was calculated from the two photon exchange diagrams by Rizov, 
Todorov and Aneva in the Quasipotential Theory for spin- 1/2 - spin- particle systems. 
The perturbation theory in the Quasipotential Theory is very similar to that of Constraint 
Theory [Eq. (2.10)], except that in the latter case off-mass shell scatternig amplitudes are 
considered and nowhere extrapolations from on-shell quantities are used. These differences 
do not, however, show up at O(a^). Furthermore, V2 is a spin independent quantity to this 
order. One finds V2P = —{V2^^y. (We also checked this result from Eq. (2.10) in the 
two-fermion case.) One then obtains for V^''^'* the following expression: 

^(1) = -(V^y + 2x2 ^(0) yio), (5^9) 

Therefore, the Hamiltonian (5.6), together with the expression (5.9) of V2^\ provides a 
gauge invariant wave equation to order c~^. Its expression in the Feynman gauge is obtained 
by setting in Eq. (5.6) U^°^ = 0. 
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The Breit Hamiltonian is obtained by using the Coulomb gauge. In the present three- 
dimensional formalism and in the cm frame, the expressions of the photon propagator in the 
Coulomb gauge and in the Landau gauge coincide in lowest order [Eq. (AlO)]. The latter 
corresponds to the choice f/*-*^-* = — ^'■^V^- By replacing U^^^ with this expression in Eqs. 
(5.6) and (5.9), one recovers the Breit Hamiltonian [^. (Notice that the term Af/'-^-'x.p 
vanishes, for AU^^^ is now a 5-function. For the same reason the term AV^'^^S.x.)^ /x."^ 
can be replaced by AV2*^°^S^/3 .) In this gauge all quadratic terms in ¥2^^ (replaced by 
its Coulombic expression (A9) ) disappear. This explains, once more, why in the Coulomb 
gauge one obtains the correct 0{a^) effects with one photon exchange diagrams only. 

The reconstitution of a covariant electromagnetic potential from the previous results is, 
however, model dependent and sensitive to implicit assumptions made about the contri- 
butions of the multiphoton exchange diagrams. A first model is provided by the Todorov 
potential [1^, introduced in Quasipotential Theory and later used in spectroscopic cal- 
culations (for electromagnetic as well as for short distance quark interactions) by Crater, 
Van Alstine and collaborators 0,^. In our notations, it corresponds to the choice (in the 
Feynman gauge) : 

1 2^^°^ 
V2 = U, = -ln(l + -^), T4 = 0, 

4 -TL 



^2^°^ = - , r = J-x^^ . (5.10) 



r 



It was shown in Ref. 0] that for large values of the coupling constant {a ~ 1/2) this 
potential leads to an instability of the vacuum state and might play a role in mechanisms of 
spontaneous chiral symmetry breaking. 



A second potential was also considered in Ref. |2^; it corresponds to the choice (in the 
Feynman gauge) : 



V2 = = ^- , T4 = , (5.11) 



1 V2' 



¥2'^^ being the same as in Eq. (5.10). This potential has the property of being regular at 
the origin and does not lead to an instability of the vacuum state. 
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C. Scalar potential 



We now turn to the case of scalar potentials and, for simplicity, shall ignore the other 
potentials. We first consider exchange of massless particles, v/^"* has then a Coulombic 
expression (A4). The term v}^^ can be calculated in perturbation theory from the two- 
particle exchange diagrams [Eq. (2.10)]. One finds : 

y« = , (5.12) 

with v/°^ given in Eq. (A4). (The corresponding calculations will be reported elsewhere.) 
As in the electromagnetic case, it is possible to modify the wave function representation by 
means of the inverse of the transformation (5.7), with U^^^ now replaced by Vi^^ /2 , and to 
reach a Coulomb-gauge like representation. It can be checked that in the resulting Hamil- 
tonian all the quadratic terms in the potential disappear. Also, the momentum dependent 
part of the effective potential then coincides with that obtained by Olsson and Miller ||25|| , 



who, with a different method took into account the retardation effects. (For a discussion 
see also Ref. PB|.) 

The foregoing results lead us to the general conclusion that, when covariant propagators 
are used, it is necessary to evaluate the leading contribution coming from the two-particle 
exchange diagrams, in order to obtain the correct 0(c^^) effects. It is only in the Coulomb 
(or Landau) gauge for the electromagnetic interaction or in a particular representation of 
the wave function for the scalar interaction that the above evaluation may be circumvented. 

As far as the spin dependent terms are concerned, Hamiltonian (5.3) agrees with the 
general expressions obtained previously in the literature . 



Scalar potentials play also an important role in the description of interquark confining 
interactions. Confining potentials, however, do not correspond to simple diagrams of pertur- 
bation theory and therefore they may display complicated or nontrivial energy dependences, 
other than those of Eq. (5.1), implied by the underlying dynamics. These energy depen- 
dences do not affect, to order 1/c^ , the spin dependent part of the Hamiltonian, but modify 
its momentum dependent part. 
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To analyze the structure of the Hamiltonian corresponding to a confining scalar potential, 
we shall adopt for the latter a more general energy dependence than that given in Eq. (5.1) 
by perturbation theory : 



1 



(0) 



^ (Pl - M) dVl 



(0) 



Pj=M 



V}°^ . (5.13) 



Here, is the nonrelativistic confining potential, while the second term in the right-hand 
side of the first equation takes into account the leading energy dependence of the relativistic 
potential Vi . As in the perturbative case, v}^^ is assumed to be a quadratic function of V^^^ 
and presumably arises from the expansion of the relativistic potential Vi in terms of v}^^ . 

The expression of the Hamiltonian is dependent on the wave function representation, 
which is modified by transformations of the type (5.7). The latter are equivalent to the 
introduction of pure gauge vector potentials. These are described by the potentials U4 and 
T4 , the leading terms of which are parametrized as : 

C/4 = — C/(°^ , T4 = — T^'^) ; (5.14) 

they satisfy the second of Eqs. (5.5) : 

Notice, however, that U^^^ and T^^^ may still depend on ratios of masses, not present in 
perturbation theory. 

As shown in subsection B, pure gauge vector potentials do not contribute to the spin 
dependent terms (to order 1/c^). 

We shall ignore in our analysis the short distance vector potential, which has the structure 
of the electromagnetic potential, discussed in subsection B, and does not affect the leading 
large distance properties of the confining potential. 

The resulting Hamiltonian takes the form (in the cm. frame) : 
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4 m\ ml mim2 ' 4mim2 ' 



2mim2 ' ^ mim2 

8 mf m2 4mim2 



(5.16) 



(The prime represents derivation with respect to and A is the Laplacian; [ , ]+ is the 
ant icommut at or . ) 

A possible comparison of this Hamiltonian with other expressions obtained from quantum 
field theoretic calculations would permit the determination of its various pieces and, with 
reasonable assumptions, might lead to a reconstitution of the covariant expression of the 
confining potential. 

In this respect, we mention the work of Barchielli, Montaldi and Prosperi ||28|, who. 



continuing the Wilson loop approach to the confining potential [^, as was developed earlier 
by Eichten and Feinberg and by Gromes [0, obtained from QCD an expression for the 



Hamiltonian to order 1/c^, including its momentum dependent part. 

The complete reconstitution of the confining potential from the comparison of the Hamil- 



tonian (5.16) with that of Ref. p8| goes beyond the scope of the present paper and is left for 



future work. We simply sketch here the qualitative features obtained from such a compari- 



son. As was already pointed out by Olsson and collaborators |^ , the Hamiltonian of Ref. 
P^l favors the interpretation of the underlying dynamics with the semiclassical flux tube 
picture. Actually, this result is not very surprising, since one of the initial hypotheses used 
in Ref. ||2^ , namely the dominance of the longitudinal color electric field, is itself motivated. 



at least partly, from that picture. 

The expression of the relativistic confining potential Vi suggested by the above compari- 
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son, turns out to be momentum dependent, typical of flux tube models. This means that in 
x-space Vi should be represented by an integral operator (in the three-dimensional variable 
x"^). However, in the classical case, the momentum dependent part is a smooth function, 
and therefore, as a first approximation, one can use for it a mean value approximation, thus 
representing Vi by a local function in . 

We also notice that the momentum dependent part of the Hamiltonian of Ref. does 
not vanish in the limit m2 oo . Therefore, in this limit, particle 1 will satisfy a Dirac 
equation with still a momentum dependent potential. 
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VI. SUMMARY AND CONCLUDING REMARKS 



We showed that the two-fermion relativistic wave equations of Constraint Theory can be 
reduced, by expressing the components of the 4x4 matrix wave function in terms of one 
of the 2x2 components, to a single equation of the Pauh-Schrodinger type, vahd for all 
sectors of quantum numbers. This equation, which is relativistic invariant, can be analyzed 
and solved with the usual tools of nonrelativistic quantum mechanics - a feature which 
considerably simplifies the treatment of two-fermion relativistic problems. 

The interaction potentials that are present in this equation belong to the general classes 
of scalar, pseudoscalar and vector interactions and are calculable in perturbation theory 
from Feynman diagrams. They are of the quasipotential type, in the sense that even in local 
approximations, which we adopted throughout this work, they may still depend on the cm. 
energy of the system. This property allows one to take into account leading contributions 
of nonlocal effects. 

In the limit when one of the masses becomes infinite, the equation reduces to the two- 
component form of the one-particle Dirac equation with external static potentials. 

The Hamiltonian of the system to order l/c? reproduces most the known theoretical 
results obtained by other methods. Furthermore, in the electromagnetic case, gauge invari- 
ance of the wave equation is checked to this order, by considering the photon propagator in 
arbitrary covariant gauges. 

A last application was devoted to the analysis of the structure of the Hamiltonian, to 
order , in the case of confining interactions. We emphasized here the role of the cm. 
energy dependence of the relativistic potential. We displayed, for confining scalar potentials, 
the general expression of the corresponding Hamiltonian, in the presence of pure gauge vector 
potentials. This expression can be used for comparison with other theoretical evaluations 
and reconstitution of the relativistic confining potential. 

The relativistic invariance of the reduced Pauli-Schrodinger type equation also allows 
one to consider ultra-relativistic limits, which play a crucial role in particle physics. 
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APPENDIX A: THE POTENTIALS FROM ONE-PARTICLE EXCHANGE 

DIAGRAMS 



In this appendix we display the relationships of the various potentials, introduced in 
Sec. II and used throughout this paper, with the propagators of the exchanged particles 
in perturbation theory. According to formula (2.10), the three-dimensional reduction of a 
propagator with mass /i results, in momentum space, from the operation : 

D{q^,l^) = D(gL = 0,/,/x) , (Al) 

or in x-space : 

D{x^,l^) = [ dxL D{x,fi) , (A2) 



D representing here the usual four-dimensional propagator. For scalar particles and in lowest 
order of perturbation theory the propagator is D{q, /i) = i/ (g^ — fi^ + ie) . 
For scalar interactions one has : 

V,{x^) = {-if ^gig2D{x^.t^). (A3) 

gi and g2 being the coupling constants of the external particles to the exchanged particle 
with mass In lowest order : 

For pseudoscalar interactions, the potential V^, has the same type of definitions as Vi . 

For QED one has : 

= -( g^,D{k, 0) - k^KF{k) ) , (A5) 
where the gauge function F is taken as a covariant function of k . Then : 



D^,{x') =-{g^,D{x\Q) + dl,diF 



X 



-( g,,D + 2gl^F + 4xJx^F ) , (A6) 
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the dot operation being defined in Eq. (2.24). 
By identification with Eq. (2.20) we find : 

T4 = 2x^2 U . (A7) 

For particles with opposite charges we have : 

V, = ^^e'D{x^,0) . (A8) 

In lowest order, V2 becomes : 

= r = a = ^. (A9) 

(Notice that the sign of the charge of the antiparticle is taken into account by the matrices 
72 of Eq. (2.20) which act on the wave function from the right. In the nonrelativistic limit, 
the matrix 72L has eigenvalue —1 .) 

For covariant gauges described by a parameter ^, we have in lowest order : 

D,Ax^) = -( g,. - f + ^ f ) D{x^, 0) . (AlO) 

The Feynman gauge is obtained with = and the Landau gauge with ,^ = 1 . 

For vector interactions resulting from the exchange of a vector particle with mass /i , the 
function F in Eq. (A6) should be replaced hj D/fi"^ . 

As a final remark, we emphasize that relationships like (A3) and (A8) between the po- 
tentials and the exchanged particle propagators are not equivalent to the instantaneous 
approximation. The perturbation theory expansion of Eq. (2.10) always reproduces for the 
potential corresponding to a one-particle exchanged diagram (including vacuum polarization 
effects) a local expression in x"^, together with the energy factor 1/{2Pl) . In the present for- 
malism, the instantaneous approximation produces a nonlocal potential due to the presence 
of the momentum operator p^"^ in it On the other hand, the instantaneous approximation 
does not provide the correct infinite mass limit, neither the correct expansion. 
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APPENDIX B: USEFUL FORMULAS FOR THE 7-MATRICES AND THE SPIN 
AND ORBITAL ANGULAR MOMENTUM OPERATORS 



In this appendix we present formulas for the 7-matrices and the spin and orbital angular 
momentum operators that are used during the reduction operations from Eqs. (2.21) to Eq. 
(3.20). We also exhibit the relationships between the components ipi , ip2 and 1^4 with ip^ 
[Eq. (3.1)]. 

The transition from Eqs. (2.21) to Eqs. (2.23) is obtained with the use of the following 
properties of the 7-matrices : 

e-^^ ill ± tL) = e±2At/4 [ (^t ^ ) 

+ i^lx^ ± ^Ix^) -L (e±2AT4 _ 1) ] , (Bl) 

where V is defined in Eq. (2.20) and V is defined as : 

V = 2Vi-V . (B2) 

Equations (Bl) are easily established by differentiating with respect to the parameter A . 

We notice that all the bilinear forms of the 7-matrices present in V [Eq. (2.20] commute 
with each other. 

Similarly one establishes the following relation giving the transformation law of under 
the action of the exponential of V : 

^-XV pT ^XV ^ pT + 2lhX [ V, + 7157251^3 

T T 

+ 7^72 ( 9ti:V2 + gl^iJ^ + ) ] xl 

+ ^ [ (sinh2Ar4) hlall^^ + lLli-x^) 
+ z(cosh2AT4 - 1) (a?33 + aJJ^) x^^ 
- 2(sinh2Ar4) 4j li.x^l2.x'^xl ] , (B3) 



where the dot operation is defined in Eq. (2.24). 
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The transition from Eq. (3.17) to Eq. (3.20) uses the following properties of the spin 
operators, defined in Eqs. (3.3)-(3.5) : 

-A + wB iw.s^W.sype^ + '^B = e'^'^^ [{W,s±W,s).p 

+ 2ih {A^2B- wB) {W^s ± W2s).x^ ] , 
-A + wB i^w,s±W,s). e^ + ^^ = 6^25 i^w,s±W,s)a , 

- 4j (1 - e'FSC) (iy,5 ± 1^25) .a; , 



+ (1 - e^C^) {W^s - W2s)a {Wis - W2s).x^ + n^P^ (1 - cosh2C) xl 



- 4(sinh2C) xl + 8C W,s.x^W,s.x^ xl ] 



(B4) 



Here A, B, C are functions of (and eventually of P^). 

The operator w can be decomposed along projectors, Wq and Wi, on spin and spin 1 
states, respectively, which also satisfy with Wi2 the following relations : 

w = Swq-Wi, Wq = ^(1 + w), wi = ■^(3-w), 

W(l-Wi2) = -{l-Wu), W(l + Wi2) = 2W + 1 - ^12 , 

Wo ^(1-^12) = 0, Wo ^(1 + ^12) = Wo, 

Wli(l-Wl2) = i(l-Wi2), Wii(l+Wi2) = ^(l+Wi2)-Wo, 



(B5) 



We also list here several useful formulas (in the cm. frame) involving the momentum, 
the spin and the orbital angular momentum operators [Eqs. (3.18)] : 
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[ p^ / ] = - ( mf x.p + 6h'f + Ah^y?r ) , (B6a) 
[p\ Sug] = -Su{ x.p + 6n^g' + Ah^^'^g" ) 

+ ^( 4:Si2 L.S - 8S1.S2 + 6h^Si2 ) g (B6b) 

(/ and g are functions of x^ and eventually of and the prime operations are defined in 
Eqs. (3.21)), 

(si ± S2).x (si ± S2).p = ^(1 ± '^12) ( ^^x.p + ihL.S ) , 

ih 

(Si ^ S2).X (Si ± S2).p = T "^12) L-(Sl - S2) , 

[ S^2 , L.(si - S2) = , [Su, [ , L.(si - S2) ] ]+ = , 
w{l + Su) L.(si - S2) = -(1 + ^12) L.(si - S2) , 
w (1- S12) L.(si - S2) = 3(1 - 5*12) L.(si - S2) , 

(l±5l2) L.(Si-S2) (l±5l2) = 0, 

(1 ± Su) L.(si - S2) (1 T Su) = 2L.(si - S2) ± [ ^12 , L.(si - S2) ] , 

[ ^12 , L.S = 2S^ - 3^2 - 3^^5i2 , 

(1 + ,^12) L.S (1 + ^12) = -2;i'(l + 5i2) , 

(1 - ^12) L.S (1 - ^12) = 4^2(1 + ^12) - 4S\ 

(1 ± >Si2) L.S (1 T ^12) = 2(L.S + S2-%^) ± [S,2, L.S], 

n x^ 

[ ^12 , [ ^12 , L.S ] = , ^12 [ ^12 , L.S ] = 2 ( L.S + S=^ - ) 

(B7) 

([ , ]+ is the anticommutator), 

ih 1 
4 si.p S2.p = h' S12 P' + ^ [ 5i2 , L.S ] x.p - — [ S12 , L.S ] L.S . 

(B8) 
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The relation expressing ■04 in terms of ■03 is : 



V'3 



_ ^-[W, + 2(1^, + 1^2) + 2/.] r 1 [ ^ L.S ] ( -z;. x.p + L.S 

I 2ai x^ 



+ ^ (1 - [ ^12 , L.S ] x.p 

- [2{Vi-Vi + Vi + U',)e-'^T, _ ^(l-e-2^4)](L.S + S^-^) 



+ \T'e 



-2Ta 



(B9) 



Notice that the terms in the curly brackets disappear in the sectors with quantum num- 
bers j — i . 

The relation expressing tjj^ in terms of ip^ can be deduced from Eq. (B9) with the 
following substitutions : ■^a <-> , ^ P'^ , <-> and a global change of sign in front 
of the curly brackets. 

Equations (3.14), which express ipi and ip2 in terms of , can be rewritten in more 
exphcit forms when potential V is given by Eq. (2.20) : 

2 



X e 



-2(t/4 + V2 + /.) I _ ,^).p _ ^ (1 _ ,-2T4) _ 



+ 



ih ,^ _2T4^ ^. 2S^ 



2x2 



(1-e- 



(4 



) (Sl - S2).X 



2S^ 



?7.' ) e 



-2T4 



(Si - S2).X 



[m 



M2 

+ §(l-e-2^)S.x 



S.p - ^(1 - e~2^4) s.x x.p 



+ 2i;i ( Vi - Vi + Vi + C/^ + Ti ) 6-2^4 s.x I V'3 , 



(BlOa) 
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i>2 



hM 



^ -2{U, + 14 + I _ s.p + -L (1 - S.X x.p 



^(l-e-2r^) S.X 



2ih ( VI - + VI + U'^ + T'^) e~2^4 S.X 

2S2 



(m^ - ml) 



+ :^ (1 - (4 



(si-S2).p 2^4) (si - S2).x x.p 



X' 



2x2 



-) (Sl - S2).X 



r)02 

2i;i ( V^/ - VI + Vl + U'^-Ti-—^U'^) (g^ _ ^^y^ 

1 1 



i>3 ■ 

(BlOb) 
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APPENDIX C: CANCELLATION OF HIGH ORDER LADDER DIAGRAMS IN 

THE INFINITE MASS LIMIT 



We show in this appendix that in the infinite mass hmit (m2 oo) the high order 
ladder and crossed ladder diagrams are cancelled by the "constraint diagrams" , so that the 
dominant contribution comes from the one-particle exchange diagrams only. 

We shall ignore the contributions of vertex corrections, which yield nonlocal terms in 
the potential, even in the one-particle problem. For simplicity, internal fermion propagators 
will be considered as free, but this assumption does not restrict the generality of the result. 
Vacuum polarization may be included in the exchanged particle propagator. For definiteness, 
we consider the case of QED (with coupling constants ei and 62 for fermions 1 and 2, 
respectively), but the derivation can be repeated with other types of interaction as well. 
Mass of the exchanged particle is irrelevant. 

In the approximation of ladder and crossed ladder diagrams, the off-mass shell (ampu- 
tated) scattering amplitude can be decomposed as : 

00 

T = E^n' (CI) 
n=l 

where is the partial amplitude corresponding to n exchanged photons. A similar decom- 
position also holds for the constrained amplitude T (2.10). 

The iteration of Eq. (2.10) yields for the potential the expression : 

00 00 
V ^ f ^{Goff ^ (C2) 

p=0 n=0 

where is that part of the potential which comes from the contributions of n exchanged 
photons. The first two terms of V are : 

Vi^f,, (C3) 
^2 = T2 + fiGoTi . (C4) 

The corresponding Feynman diagrams are represented in Fig. 1. While T2 contains the two 
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usual Feynman diagrams of two exchanged photons, the last term in Eq. (C4) yields the 
"constraint diagram" (third diagram in Fig. lb). 

The expression of Vi was given in Appendix A (notice that to this order the potentials 
V and V [Eq. (2.17)] coincide); in the limit 1712 oo , Vi behaves as 0(l/m2). We now 
concentrate on the behavior of V2. 

In T, the external momenta are submitted to the constraint (2.11) and, therefore, satisfy 
the relations : 

Pl+P2= P'l +P2 = P , 

2 2 /2 /2 2 2 / / ^r'tx\ 

Pi - P2 = P 1 - P 2 = ^^-^2 , PiL=PlL , P2L= P2L ■ (C5) 

In the limit m2 — > 00, one has : 

1 1 

P20 = m2 + 0{ — ), Pq = m2+pio + 0{ — ), 

1712 ^2 
IPI 

Pio-p'io = P20-P20 = 0{ — ). (C6) 

1712 

The expression of T2 is : 

i f d^ki 



T2 = {-i)\eie2y (^) J ^Mii^i(^i) Df^2U2{Pi - p'l - h 



XTlMi-^llPl - ^l)7l/.2 [ 72i.i^2(-(P2 + ki)h2u2 + 12V2S2{-{P'2 " ki))l2u^ ] , (C7) 

where the photon propagator is considered in an arbitrary covariant gauge [Eq. (A5)]. 



In the limit ra2 ~^ 00 one has 17,18 



[ l2u^S2{-{p2 + ki)h2u2 + 12U2S2{-{P2 - fcl))72z.i ] 

1 % 

= [ 721^1 7 —r-^ -—721^2 + 721^2 Tl ^-^ T— 72!/! ] 

-72- (P2 + h) -m2 + ie -12-{P2 -h) -m2 + te 

= (-l)'5.io5.20 27r6{ho) + O(^) , (C8) 

where the eigenvalue —1 has been used for the matrix 720 (antifermion). (Because of the 
factor 6{kiQ) and Eqs. (C6), the gauge part of the photon propagator disappears in Eq. 
(C7) at leading order.) 
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On the other hand, in the iteration term TiGqTi [Eq. (C4)] the integration concerning 
Go is three-dimensional after the use of the constraint (2.11) with the fermion momenta 
(pi + ki) and (p2 — ^i)- This yields the condition kn = . Hence, one obtains : 



^lliilSi{pi - /Cl)7l^2 l2u^S2{-{p2 + ki))l2u2 Hq{P2 + ^l) • (C9) 

Notice that because of the constraint (2.11) Hq{pi — ki) = iJo(P2 + ^i) [Eq. (2.13)]. One 
also has : 

S2{-{p2 + ki)) Ho{p2 + ki) = i{ -72.(^2 + ki)+m2) . (CIO) 

In the limit m2 oo , this term, multiplied by the matrices 721/^ and 72,^3 , behaves as 
(— l)^5,^jO<^j^2o(^2m2) and the factor 2m2 cancels one of the factors 2Pl of the denominator 
of the right-hand side of Eq. (C9). Furthermore, in the same limit S{kiL) — > 6{kio) and, 
because of the additional factor i"^, the iteration term TiGqTi cancels, at leading order, the 
amplitude term T2 in Eq. (C4). Hence, V2 behaves, when m2 — > 00 , as 0{l/m'^ . 

The above calculations can be repeated at higher orders of perurbation theory. For the 
diagrams with n exchanged photons, Vn is the sum of the partial amplitude T„ and the 
iteration terms coming from VGqT ; it can be written in the form : 

Ki = + X/ X/ ^riGoTfjCo ■ ■ ■ T^pGoT^p+l , (Cll) 
p=l riH hrp+i=n 

where, in the generic term of the sum, the constraint factor Gq appears p times. A typical 
diagram, where p = 2, is shown in Fig. 2. 

Putting aside all factors which do not play an essential role, an amplitude , that 
contains all ladder and crossed ladder type Feynman diagrams with r exchanged photons. 



behaves in the limit m2 00 as p]7|JT 



T, ~ {^-^f^{^-e,e2r (-^) (27r)'-i 5{k,,) ■ ■ ■ 5{K-^,) , (C12) 

while the occurrence of a constraint factor Gq just after provides the contribution : 
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Go ~ («2m2) 2nS{kro) ■ (C13) 

It can be checked , using these results, that a term of the iteration sum in Eq. (Cll), 
containing p factors Go , behaves as : 

Tr^GoTr^Go- ■ ■TrpGoTrp_^_^ ~ (-1)^T„. (C14) 

Furthermore, the p factors Go may appear in such a term in ("p~^) different configurations. 
The behavior of Vn is then : 

_ n-l 

hm m2Vn = hm maT, [ 1 + ^(-1)^ (V^) ] = 0, n>2. (C15) 

7712— >oo m2^oo ^ — ' \ -f^ / 

p=l 

Therefore, Vn behaves as 0(l/m2) and the dominant contribution to V [Eq. (C2)] comes 
solely from the one-particle exchange diagram, which behaves as 0(l/m2) . A similar result 
is also obtained in other approaches, based on a three-dimensional formulation of the two- 
body bound state problem |1I3| , |22[| . 
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FIGURES 

FIG. 1. The diagrams contributing to Fi , (a), and 1^2 , (b) [Eqs. (C2)-(C4)]. The third diagram 
in (b) is the "constraint diagram" . 

FIG. 2. A typical "constraint diagram" of high order [Eq. (Oil)], where the constraint factor 
Go [Eq. (2.12)] appears twice (boxes with a cross). The boxes without a cross represent ordinary 
ladder and crossed ladder type Feynman diagrams (contributing to the partial amplitudes Tj. in 
Eq. (Cll)). 



49 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9403232vl 



